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. Abstract 

\D 

I/*} I We study various dynamical properties (winding angles, areas) of a set of har- 

monically bound Brownian particles (monomers), one endpoint of this chain being 
. kept fixed at the origin 0. In particular, we show that, for long times t, the areas 

<^ \ {Ai} enclosed by the monomers scale like i 1//2 , with correlated gaussian distribu- 

tions. This is at variance with the winding angles {9i} around fixed points that 
scale like t and are distributed according to independent Cauchy laws. 

a 

O 

In this paper, we will study the planar motion of a chain of n harmonically bound 
brownian particles. This model is usually refered in the litterature as the Rouse chain 
^ [1] and has shown to be historically very important in polymer science [2, 3]. We will 

consider such a chain attached at the origin and examine some of its properties from 
the Brownian motion viewpoint. 

Representing a given configuration of the chain by a complex n-vector z (the com- 
ponents the complex coordinates of the particles), we consider the 
set of all the closed trajectories of length t, i.e. z(t) = z(0), and this for all the starting 
configurations z(0). Practically, we will not weight the starting configurations with any 
thermodynamical factor. We are aware that this approach is quite different from the one 
taken in polymer physics [4] where, at t = 0, the chain is supposed to be in equilibrium 
with the environment at some finite temperature T. 
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Aj and Oj being the area enclosed by the j th particle and its winding angle around 0, 
our goal is to compute the joint probability distributions P({Ai}) and P({8i}) for such 
trajectories. In order to make comparisons, we now recall some of the results concerning 
the planar Brownian motion. 



We first quote the area and winding angle distributions, respectively P(A) (Levy's 
law [5]) and P(0) (Spitzer's law [6]) for a particle allowed to wander everywhere in the 
plane: 



P(0) = -4 9 (2) 

(the last one holds, in the limit t — > oo, for open curves, the final point being left 
unspecified). 

Those two laws were obtained more than 40 years ago and since that time many 
refinements have been brought. For instance, in [7], the authors pointed out the im- 
portance of the small windings occuring when the particle is close to 0. Excluding an 
arbitrary small zone around 0, they showed that the variance (6 2 ) becomes finite in 
contrast with the Spitzer's result, eq.(2). 

On the other hand, for Brownian motion on bounded domains [8, 9], the scaling 
variables in the limit t — > oo, become, resp., A/y/i and 9/t with still an infinite variance 
(0 2 ). We close here this brief recall and start our chain study with the following set of 
coupled Langevin equations : 



zi = k (z 2 - 2zi) + r]i 

zi = k (z l+1 + - 2zi) + rji , 2 < I < n - 1 (3) 
z n = k (z n —i — Z n ) + T) n 

where k is the spring constant and r\ m (= r\ mx + ii] my ) a gaussian white noise: 



(Vm(t)) = 

(rjm(t)Vm'(t')) = 2 5 mm ,5(t-t') (4) 
Introducing the complex n-vector r], eq. (3) can be written in a matrix form: 



z = — kM.z + 7] 
2 



(5) 



where M is the tridiagonal (n x n) matrix: 
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with an inverse given by: 
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The eigenvalues of M are: 

Uj = 2 1 - cos - , 1< j < n 

With the matrix u = diag(cjj), we can write: 



u = R^MR (7) 

z = KZ (8) 

where R is an orthogonal matrix and the components of Z are the normal coordinates. 

Let us call V(z, z , t) the probability for the chain to go from z at t = to z at time 
t. V satisfies a Fokker-Planck equation [10]: 

d t V = [ t d z kMz+ t d,kMz + 2 t d,d z )v (9) 

where d z (resp. d z ) is a n- vector of components d Zi (resp. d Zi ) and t d z (resp. t d z ) is 
the transpose of d z (resp. d z ). The solution can be written in terms of a path integral 
(VzVz = U?=i'Dz i Vz i y. 

V(z,z ,t) = det (e' feM ) / PzPzexp ( — / '(I + feMz)(i + fcM^) dr 1(10) 
v / j 2 (o)=2 V 2 jo / 

= F(z,Zo,t).G(z,zo,t) 

with 
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G(z,z ,t) 



F(z,z ,t) 



H 




-tH 



^o)(H) 



(12) 



(13) 



The matrices S and C appearing in (12) are defined as: 



S = kM (sinh(tfcM)) 



-i 



C = kM coth(tA;M) 



(14) 



In fact, V, eq. (10), can be easily deduced from the gaussian distribution of r\ (use 
(5); det(e* fcM ) is simply the functional Jacobian for the change of variable rj — > z [11]). 

(12) is a generalization of the harmonic oscillator propagator [12]. It is obtained 
by using the normal coordinates. Furthermore, as can be easily checked, V is properly 
normalized: / dzdz V(z,z ,t) = 1. 

Remark that an effective measure can be built for a distinguished monomer of the 
chain [4]: this can be done by integrating the Wiener measure (10) over all the paths of 
the other monomers. The result is a complicated expression that contains, in particular, 
a non local part (in time) exhibiting the non-Markovian character of the process for this 
monomer. Nevertheless, we will show, in the sequel, that, despite this complication, we 
can compute some joint laws for several monomers (and a fortiori for one monomer). 

So, let us turn to the computation of the area distribution P({Ai}) for closed tra- 
jectories. Inserting the constraint 



n particles submitted to uniform magnetic fields orthogonal to the motion plane (in 
addition to the harmonic interactions). Remark that, in principle, the magnetic fields 
are not the same for all the particles. 

Introducing the (n x n) diagonal matrix B (B^ = BiSij), we obtain 




(15) 



in the Wiener measure 





(16) 
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Z B (t) = Tre- tHB 
H B = H + V 

V = I PzB9 z +^-B9,) + l t zB 2 z (17) 

In general, the matrices B and M do not commute and it is a difficult task to get 
the partition function Z B (t). On the other hand, the distribution of the total area 
A = Yh=i M is obtained by taking Bj = B for all j. In this case, B and M commute. 
Using normal coordinates and known results about the partition function of the "2D 
harmonic oscillator + uniform magnetic field" problem [13], we get the characteristic 
function of A (I n is the (n x n) unit matrix): 



M = TT ( COShjtkVj) - 1 , 

Zo(t) j-Ji Vcosh(i v /(^) 2 + (f ) 2 ) - COSh(if) ' 

det (cosh(tA;M) - I n ) 



det (cosh(t v /(A;M) 2 + (f ) 2 ) - cosh(tf)) 
(the u)jS are defined in (6)). By Fourier transformation, the above j th factor gives: 



(19) 



«,■ 



4— sinh 

7T 
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m,r=0 

ajt(r + 1/2) +iajA 
aijt(m + 1/2) — iajA 
k ujj 




m'Jjr 



j in 



'jm IjrJ 



20) 



where the K±i are modified Bessel functions. 

Thus, in the general case, P(A) can be obtained by a n-convolution product of the 
Pu/s- However, we are afraid that the final result could be awkward! Nevertheless, if 
we consider the limit t — > oo, (18) leads to: 



z B (t) 
z (t) 



exp 



~8k 



E 



cxp 



tB 2 n(n + 1)' 
16k 



(21) 



Then, Fourier transformation shows that, in the large t limit, A is gaussian and scales 
like y/i. Such a scaling is expected for all the areas A { . This is what we will demonstrate 
by perturbation theory. When t — > oo, we have 



Z B (t) ~ e~ tE °W (22) 

where E'o(B) is the ground state energy. Moreover, due to the large oscillations of the 
factor e lBj j in (16) when Aj — > oo, only small values of Bj will contribute. So, it is 
enough to compute Eq(B) at lowest order in B. We will use the normal coordinates Z{. 
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The eigenstates of H are given by [14] 



*w,K>({^}) = n 



i y\ A n i + K'l 

n 

E {mj },{ nj } = E( 2n i + Kl + !) 



-/ m > Q i {ujj \Zj\ 2 ) ^ 1 /2 L^ 1 \Zj\ 2 )e-* u * 1 ^ 1 



(23) 
(24) 



where L^ 1 is a Laguerre polynomial and the ground state is ^{o},{o}- The pertur- 



bation V, (17), writes: 

V = - (-*ZR _1 BRd z +* ZR _1 BR^) + -*ZR _1 B 2 RZ 



(25) 



At first order in V, we get: 



a^ 1} (b) = J *; 0}>{0} ™ 



1 1 n 

{0},{0} = ^Tr(B 2 M- 1 ) = -E^^ 



2 

S/C m=l 



(26) 



Quadratic terms in B will also be produced at second order in V. The no n- vanishing 
contributions will come out from the transitions from the ground-state to the states 
{rrij = ±1, mi = =1=1, — if i ^ j, I}, {n k = 0}. The computation gives: 



A£j 2) (B) 



1 



<±n _i_ <£± 



, „ , ^ ' B rn B m i ^ ' R m iR m jR m iiR m rj 

lOft a _l, \ UJl + UJj 



m,m'=l 

1 n 



3+1 
1 n 



with: 



m=l m,m'=l 



i n p2 p2 
1 \ ~> -n-m/ -"-m'/ 



4 JO L V /m,m'J 8 

So, to lowest order in B, we get: 



Z B (t) / t " 

ZoW~ eXP r2l ^ BmD m , m ,B r . 

U V / \ m,m'=l 

As can be easily checked, (21) is recovered if we set B m = B, Mm. 
With (16), we arrive at the probability distribution: 



(27) 
(28) 



(29) 



(30) 
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™ = [It) vab ^ (-1 J_ ^-u,^) < 31 > 

Thus, we observe that the areas Ai are correlated gaussian variables and that they 
scale like t 1//2 as expected. For the special case n = 2, we have: 



P(A 1 ,A 2 ) = 




The width of A 2 is larger than the one of A^. this is related to the fact that the second 
particle is, in average, farther from than the first one. So, it sweeps larger areas. 



Now, going to the winding angles around 0, we proceed as before and insert the 
constraint 

in the Wiener measure. We are now faced to the problem of n harmonically bound 
particles submitted to the magnetic fields of point-like vortices located at the origin. 
The corresponding hamiltonian is: 



if A = H Q + W (34) 

n , -1 is n \ 2 

w = 5> ( ^_ aJ+E^V (35) 

and the distribution P({^}) is given by: 

*<w>=/(n£^)(fi) («) 

Studying the limit t — > oo, we cannot develop directly as before a perturbation 
theory with W: this is because of the last term in W that leads to a singular pertur- 
bation [9]. Due to this term, all the eigenfunctions of H\ must vanish in at least as 
n™=i(^i)' Ai '/ 2 (= U) . So we redefine those eigenfunctions [9]: 

= UV (37) 

The new hamiltonian acting on \l/ is: 



7 



H = H + W (38) 

~ n / 1 1 \ 

W = E (( A * ~ \ X i\)- d zi ~ (K + \K\)-dzi) (39) 



That time, we can compute AE (X) perturbatively and it will appear that only first 
order is necessary. Integrals of the form 

J e-^ kMz jd Zi e-^ kMz dzdz (40) 

are involved. Integrating by parts and using d Zi = ir5(zi), we get, after some 
algebra: 

n I \ I n I \ I 11 

^W = *E 7 i^ ! r = *E 1 T i = E^I ( 41 ) 
j=i i iv± )jj j=i j j=i 

So, for the winding angle distribution, we obtain: 



= n 



(42) 
(43) 



At large times, the winding angles are uncorrelated, they scale like t and are dis- 
tributed according to Cauchy laws. The variance (9^) is infinite: this is, of course, due 
to the "small windings" occuring in the vicinity of the origin as will be seen explicitly 
at the end of this paper. 

Moreover, we observe that 6j scales like /Xj, i.e. like 1/j. This is reasonable because, 
when j increases, the considered particle is, in average, farther from and, consequently, 
its winding angle must decrease. What is somewhat unexpected is such a simple depen- 
dence of 0j on j. 

We also addressed the problem of winding angles around n different points of complex 
coordinates bi, I — 1, . . . n. 

O'j being the angle wound by the particle j around the point bj, we obtained for the 
set of variables {6^} the same joint law as (43) except for the change of fij into ji'y 



(44) 



Owing to the rotationnal symmetry breaking when bj ^ 0, the winding angles 0j 
are statistically reduced by the factor e - ^'^' 2 . Nevertheless, even for large |6j|'s, the 
variance {(9'j) 2 ) is infinite. 

Setting all the b/s to zero, we recover (43). This is what we will consider now and 
assume that we count the winding angles 0j only when \zj\ > r$ (i.e. the so-called "big 
windings" [7]). Still when t — > oo, the perturbation W, eq.(35), can now be used because 
Xj = when \zj\ < r . At first order in W, the linear contributions in Xj will cancel. In 
the limit of a small, but finite r , we get, for the remaining contribution: 

n \2 

AE^\X) ~ A;|lnr |^^ (45) 
j=i J 

The quadratic contributions in the X/s coming out from the second order in W will 
be finite (thus subleading) when r — > + . Finally, we get for the big winding angles 
asymptotic distribution: 

P(R-}) = fl \l A ,n , exp (-—1 rei] (46) 

j-Ji\j 4:7rtk\\nr \ F \ 4^|lnr |V v ; 

In this limit, the variables 9j are uncorrected: the correlations get smaller and smaller 
when r decreases. They are now gaussian and scale like \Jt\ \nr Q \/j. Their variance 
grows to infinity when r goes to 0, showing the increasing contribution of the small 
windings around 0. 

To summarize, we have computed explicitly the asymptotic joint laws for the areas 
(that scale like \/t) and for the winding angles (that scale like t when no critical region 
is excluded). The scaling variables we have got compare well with those involved in the 
Brownian motion on finite domains: this is not so surprising since the chain is bound to 
a fixed point. 

Moreover, we have shown that physical interactions between particles (harmonic 
interactions here) can lead to statistical correlations (case of the areas) or not (case of 
the winding angles): it depends on the quantity we consider. 

In a forthcoming paper [15], we will study the statistical properties of the free Rouse 
chain. We will especially show that the areas and winding angles distributions are very 
different from those presented in this work. This is essentially due to the translation 
invariance that holds when the chain is free. 



One of us (O. B.) acknowledges Dr. G. Oshanin for drawing his attention to this 
problem. 



9 



References 

[1] RE. Rouse, J. Chem. Phys. 21:1273 (1953). 

[2] P.G. de Gennes, Scaling Concepts in Polymer Physics (Ithaca, 1979). 

[3] M. Doi and S.F. Edwards, The Theory of Polymer Dynamics (Oxford, 1986); A.Y. 
Grosberg and A.R. Khokhlov, Statistical Physics of Macromolecules (New- York, 
1994, translated from russian). 

[4] S.F. Burlatskii and G. Oshanin, Theor. and Math. Phys. (3) 75:659 (1988, trans- 
lated from russian). 

[5] P. Levy, Processus Stochastiques et Mouvement Brownien (Paris, 1948). 

[6] F. Spitzer, Trans. Am. Math. Soc. 87:187 (1958). 

[7] P. Messulam and M. Yor, J. Lond. Math. Soc. (2) 26:348 (1982); J.W. Pitman and 
M. Yor, Ann. Prob. 14:733 (1986); 17:965 (1989); J. Rudnick and Y. Hu, J. Phys. 
A 20:4421 (1987). 

[8] J. Desbois and A. Comtet, J. Phys. A 25:3097 (1992). 

[9] A. Comtet, J. Desbois and C. Monthus, J. Phys. A 26:5637 (1993). 

[10] N.G. van Kampen, Stochastic Processes in Physics and Chemistry (North-Holland, 
Amsterdam, 1981). 

[11] F. Riesz and B. Sz-Nagy, Functional Analysis (Frederick Ungar Publ. Co., 1955). 

[12] R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals (New- York, 
1965). " 

[13] A. de Veigy and S. Ouvry, Nucl. Phys. B 388:715 (1992). 

[14] L. Landau and E.M. Lifshitz, Quantum Mechanics (London, 1958, translated from 
russian). 

[15] O. Benichou and J. Desbois, in preparation. 

e-mail: 

benichou@lptl.jussieu.fr 
desbois@ipno.in2p3.fr 



10 



